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Abstract. Quasiclastic neutron scattering is a powerful tool for the study of 
non-periodic motions in condensed matter as a detailed line shape analysis can 
give information about the geometry and rate of the scatterers' displacements. 
Unfortunately, there are also a number of artifacts which can masquerade as signatures 
of motions and can therefore lead to erroneous results. Their influence on the 
evaluation of the motions of the phospholipid dimyristoylphosphatidylcholine (DMPC) 
is discussed. On a 60 ps time scale, the long-range motion of the molecules has a flow- 
like character with similar velocities above and below the main phase transition. It is 
proposed that the concepts of dynamical heterogeneities and "floppy modes" developed 
in glass physics provide a framework to explain the observed behaviour. 
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Figure 1. Width of the 
narrow Gaussian line of a three- 
component fit to the spectra of 
the fully hydrated phospholipid 
dimyristoylphosphatidylcholine 
(DMPC), very similar to the 
results published before p] 
where it was speculated that 
the levelling off of the widths at 
low temperatures and small Q 
could either be a signature of 
confinement of the molecules or 
an artifact. The distance between 
two data points in the spectra, 
SE = lOueV, is indicated in the 
figure. The energy resolution 
is about 60ueV, corresponding 
to an observation time of about 
60 ps. The solid lines are fits to 
the filled points in the medium-Q 
region. 

1. Introduction 

Recently, a quasielastic neutron scat- 
tering (QENS) experiment on the fully 
hydrated phospholipid dimyristoylphos- 
phatidylcholine (DMPC) was published pp. 
It corroborated collective flow-like motions 
as a novel view of the long-range motion 
of phospholipid molecules that had been 
put forward by molecular dynamics simula- 
tions [2j . A fit to the data showed that the 



component which was assigned to the long- 
range flow-like motion of the molecules lev- 
elled off at temperatures below the main 
phase transition (24°C) and low values of Q 
as shown in figure [TJ This was discussed to 
be a sign of a confinement of the molecules 
during the 60 ps observation time in a cage 
of their neighbours - or an artifact [U El S] . 

The macroscopic diffusion of phospho- 
lipid molecules is known to decrease by 
about one order of magnitude when the 
membrane is cooled below the main phase 
transition [5]. In contrast, the high-Q part 
of the broadening caused by the long-range 
motion shown in figure [TJ decreases much 
less. The confinement length that was ex- 
tracted from the position of the deviation 
of the linear behaviour at low Q to about 
9 A coincides quite well with the molec- 
ular repeating distance in the membrane. 
This made a flow-like motion of the sin- 
gle molecules until they hit the cage of the 
neighbouring molecules seem reasonable. 

However, it was speculated already 
there that this might be an artifact as 
the broadening is very small compared 
to the instrumental resolution. A caging 
effect should also not only cause a 
plateau of the quasielastic broadening but 
simultaneously an emerging elastic S(oj) 
component which was not observed. On a 
physical level, it is not clear why the flow- 
like, collective motions of the molecules 
together with their caging neighbours 
above the phase transition should persist at 
low temperatures for the single molecules 
within the cages. 

The present contribution aims to 
clarify the picture of the motions of 
phospholipid molecules on a 60 ps time 
scale below the main phase transition. The 
two likely artifacts that could simulate the 
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observed plateau of the broadening at low 
Q are studied: multiple scattering and fit 
artifacts. 

Only two years after the collection of 
the first neutron diffraction patterns ever in 
1946 [6], multiple scattering was identified 
as the source of a spurious background 
contribution [7]. A few years later, 
efforts to treat multiple scattering (semi- 
) analytically were started [81 191 [TT71 ITT] 
and later complemented by Monte Carlo 

simulations [i2i H3i da H5i dsi iizi nsi us]. 

However, still by far most quasielastic 
neutron scattering experiments study sam- 
ples with very high transmission in order to 
be able to neglect rather than correct mul- 
tiple scattering. The transmission needed 
to justify this neglect is not agreed on in 
the literature - some authors feel safe with 
samples with a transmission of 99% [20J, 
most adjust their samples to achieve a 
transmission of about 90% j2B [22] but oth- 
ers still feel comfortable with samples with 
a transmission as low as 70% [23J. 

It is also often unclear what to expect 
from uncorrected multiple scattering. An 
artificial broadening due to multiple scat- 
tering was for example assumed to yield 
a diffusion coefficient that was higher [23] 
than the generally accepted value [23 [26] . 
Others suggested that this broadening re- 
sults in a Q dependence of the line width 
as shown in figure [1] which is normally in- 
terpreted as confined diffusion [4J. Offsets 
of the line width and the elastic incoherent 
structure factor (EISF) are also often ob- 
served [TTJ [23 [211 [TJ some of which were 
shown to be due to multiple scattering. 

Whereas the neutron scattering com- 
munity is well aware of the problem of mul- 
tiple scattering, fitting artifacts are widely 
ignored which can lead to erroneous results: 



one crucial step is the numerical convolu- 
tion of the theoretical fit function with the 
measured instrumental resolution [21] . The 
instrumental resolution gets more impor- 
tant as neutron spectroscopy is employed 
for the study of very slow motions because 
the line broadening caused by the sample 
gets small. 

It is even sometimes said that the 
instrumental resolution should be seen as 
a lower limit for the detectable broadening. 
However, while this is a valid rule of thumb 
for localised motions, many measurements 
of long-range motions which cause a 
broadening of the whole line are done below 
this limit: from neutron scattering studies 
of diffusion in water [30] over proteins [2H] 
to metallic liquids [3TJ or X-ray scattering 
studies of glass formers [32] . 

In the following, a short repetition of 
the basic assumptions of neutron scattering 
data evaluation is presented, highlighting 
the points where multiple scattering and 
fit artifacts come into play. After a 
description of the studied systems, the 
effects of multiple scattering and the point 
density used for the fit are described 
generally enough to be useful also for 
other experiments. The consequences for 
the evaluation of long-range and localized 
motions are presented and with this 
knowledge, a coherent picture for the 
phospholipid dynamics is proposed. 

2. Concepts of QENS data 
evaluation 

The evaluation of QENS data relies vitally 
on the fact that the recorded number 
of neutrons function of solid angle 

and energy transfer, the double differential 
cross section, can be transformed into the 
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scattering functions S^ inc (Q,uj) via [33J 
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In this relation, it is assumed that a 
monochromatic beam of neutrons was 
scattered at most once in the sample and 
that the scattering lengths of the atoms 
do not correlate with their positions. It 
is often further simplified by neglecting 
the coherent scattering contributions of 
hydrogen-rich samples and using only the 
absolute value of the scattering vector, Q, 
for powder samples. 

To account for the non-perfect mono- 
chromaticity of the incident neutron beam 
and the finite energy resolution of the 
spectrometer, the experimentally obtained 
scattering function S^(Q,u) is under- 
stood as the convolution of the instrumen- 
tal resolution function R(Q,u) with the 
theoretical scattering function S^°(Q,u): 

S^(Q,u) = S^°(Q,co)®R(Q,u) . (2) 

For the analysis of QENS data, this 
convolution is carried out only in the 
energy-direction. As the real resolution of 
the spectrometer does not strictly follow 
an analytical function, the integral form of 
this convolution 



R(Q,u-u') du/ 



(3) 



is mostly replaced by the discrete equiva- 
lent 



N 



Sine (Q> U n) — Sfn^jQ, 



OJr. 



m=l 



R(Q,u n - u m ) ■ 5u 



(4) 



where N is the total number of points of 
a spectrum in the energy direction and Sco 
the width of a bin. 

The relevant motions of the scatterers 
are then assumed to be independent of each 
other so that the scattering function can 
be approximated composition of e. g. 
translational, rotational, and vibrational 
contributions [36] , 

sST(0,w) = fl£r(Q,w)® 

Si:l(Q,u)®S^(Q,u) . (5) 

The amount of physics and approxi- 
mations contained in these equations is so 
overwhelming that it seems safe to postu- 
late that by far most neutron scatterers (in- 
cluding the authors) settle for using these 
approximations after following the calcula- 
tions a few times. 

In this contribution, we discuss the 
influence of two of these steps: (a) The 
assumption that all detected neutrons 
were scattered only once in the sample 
and (b) the correction of the limited 
instrumental resolution by the discrete 
convolution of the theoretical scattering 
function with an experimentally obtained 
instrumental resolution function. 

3. Experimental 

3.1. Materials 

A vanadium hollow cylinder standard with 
an outer radius 11.75 mm and wall thick- 
ness 0.6 mm was measured at 20 °C. The 
transmission was calculated by numerical 
integration of the path length of the neu- 
trons through the scatterer in transmis- 
sion geometry and assuming an exponen- 
tial decay of the neutron intensity in the 
medium 1371 to 93%. 
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Figure 2. Molecular structure of 
DMPC. 

For all samples, hollow cylindrical 
aluminium sample containers j38j [39] with 
an inner radius of the outer aluminium 
layer of 11.25 mm and varying radii of the 
inner aluminium layer were used. The 
choice of the insert determined the gap size 
and therefore the sample layer thickness. 

The phospholipid dimyristoylphos- 
phatidylcholine (DMPC, cf. figure [2]) 
was obtained from Lipoid GmbH, Lud- 
wigshafen (Germany) and used as received. 
It was hydrated at 40 °C via the gas phase 
with D 2 99.90% from Euriso-Top, Gif- 
sur-Yvette (France) for at least 48 h until 
a clear phase formed. Subsequently, fur- 
ther D 2 was added up to 50 weight % 
to ensure full hydration during the whole 
experiment. With a gap size of 0.2 mm, 
a calculated transmission of 85% was ob- 
tained. The measurements were performed 
at 10 °C, 20 °C, 30 °C, and 40 °C. Also pure 
D 2 was measured at these temperatures. 

Pure H 2 0, already previously used 
for multiple scattering studies [40j |4"T] . 
was taken from a Barnstead EASYpure II 
ultrapure water system with a resistivity 
of p > 18.2Mfi-cm and filled into 
the sample containers with gap sizes of 
0.1mm, 0.2 mm, 0.3 mm, 0.5 mm, 1.0 mm, 
2.0 mm, and 3.0 mm. These samples have 
calculated transmissions of 85%, 72%, 62%, 
46%, 23%, 6.4%, and 2.2%, respectively - 
the mean free path length of a neutron in 



H 2 is about 2 mm. The measurements 
were performed at 20 °C. 

3.2. Measurements and data treatment 

All experiments were performed at the 
cold neutron time-of-flight spectrometer 
TOFTOF H2J S3] at the Forschungs- 
Neutronenquelle Heinz Maier-Leibnitz (FRM II), 
Garching bei Miinchen, Germany. The 
spectrometer is equipped with detectors in 
an angular range from 7.5° to 140°. 

An incident neutron wavelength of 
6 A was used, hereby giving for quasi- 
elastic scattering access to Q values 
from 0.3 A -1 to 1.8 A -1 . The maximal 
momentum transfer of elastically scattered 
neutrons at this wavelength is h-2.1 A" 1 for 
backscattered neutrons. An instrumental 
resolution of about 60 ueV (full width 
at half maximum) was reached with a 
chopper rotation speed of 12000 rpm which 
required a frame overlap ratio of 4. The 
corresponding observation time is about 

60ps [g. 

Today's normal way of histogramming 
the data online during the experiment 
adds a restriction to the rebinning of the 
data. The online histogramming was done 
into 1024 time channels, resulting in a 
width of about 7.5 ueV of the bins at zero 
energy transfer. These data points were 
then rebinned into bins with a constant 
step in energy of 6E=10\ieV during 
data reduction into slices of constant Q. 
Additionally, the phospholipid was also 
measured with an online histogramming 
into 4096 time channels, resulting in a 
width of about 1.9 peV at zero energy 
transfer and a constant step of <5E'=2ueV 
after rebinning. 

The data were normalized to the scat- 



Slow short-time phospholipid motions, multiple scattering, and fitting artifacts 



6 



tering of the vanadium standard to correct 
for different detector efficiencies. An empty 
can measurement was subtracted from all 
measurements and the D 2 spectra were 
subtracted pro-rata from the ones of the 
phospholipid at the corresponding temper- 
atures. The obtained spectra were not nor- 
malised to the amount of sample in the 
beam or the total scattered intensity. 

Diffraction patterns da/dQ(26) were 
obtained by integrating the double differ- 
ential cross section over all energies. Spec- 
tra were obtained by regrouping the data 
into slices of constant Q [37]. In energy, 
the spectra were evaluated in the region of 
AE = hoj from —1 meV to +1 meV where 
positive values for the energy transfer de- 
note an energy gain of the neutron during 
the scattering process(es). 

The fits of the theoretically calculated 
scattering functions S^^°(Q,u) to the 
measured ones S^(Q, u>) were numerically 
convolved with the instrumental resolution 
R(Q, uj) as determined by the measurement 
of the vanadium hollow cylinder standard. 

In order to enhance the visibility 
of the inelastic part of the spectrum, a 
representation as the imaginary part of the 
dynamic susceptibility is also shown, the 
transformation was done using 

X "(Q,u) = S e £(Q,uj)/n(uj) (6) 

where the scattering functions are divided 
by the Bose occupation factor 

n(co) = l/(exp[hw/(k B T)}-l) . (7) 



an overall scaling factor a was included into 
the fitting function. 

3.3. Fitting models 

The data of the phospholipid were fitted 
with a convolution of a Gaussian for 
the long-range flow-like motion with two 
localized diffusive contributions, each of 
the form 

S^°(Q,u) = MQ) ■ 8(f») + 

(l-A (Q))-L(T hroad (Q),co) (8) 

where L denotes a Lorentzian 

r(Q)/7r 



L(T(Q),u) 



(9) 



r(Q) 2 + w 2 ' 

The different localized contributions were 
assigned to the motions of the atoms in 
the head and tail groups, respectively []]. 
The Gaussian then captures only the long- 
range motion of the whole molecule. The 
Gaussian line shape with a line width of 

a(Q) =v-Q (10) 

results from the assumption of flow motions 
with Maxwell-Boltzmann distributed flow 
velocities around the most likely velocity 
v P. 

The spectra of H2O were fitted with 
a common function in quasielastic neutron 
scattering, the sum of two Lorentzians 
which is motivated [44J by the convolution 
of a long-range diffusive motion 

S^°(Q,co) = L(T„(Q),u) , (11) 



with a localized motion With a scaling 
Intermediate scattering functions /^(Qp^factor a accounting for the scattering 
were obtained by cosine Fourier transform 
of the energy spectra and subsequent divi- 



sion of the sample data by the instrumental 
resolution as determined by the vanadium 
measurement. Again, the data were not 
normalised to I^(Q,t = 0) = 1. Instead, 



power of the sample, this yields 
S£f° (Q, w) = a ■ L(T„(Q),u) <g> 

[MQ) ■ + 

(l-A (Q))-L(T hroad (Q),u)] 



MQ) ■ £(W)r 



Slow short-time phospholipid motions, multiple scattering, and fitting artifacts 



7 



a ■ (i - a (Q)) ■ L(r(g) narrow + 
r(g) br0 ad,w)) . (12) 

A specific model for the localised motion, 
giving an expression for Aq(Q) as used in 
previous works [30J was avoided because of 
the distortion by multiple scattering [36j 

nan]. 

For the evaluation of the intermediate 
scattering function, another fit function 
will be used which is also very common for 
the evaluation of QENS data: a stretched 
exponential, 



J i 1 n h c eo (Q,t) = a-exp -(t/r 



(13) 



The alternative, a sum of two exponential 
decays, gives very similar results to the 
ones obtained from the fit of the scattering 
function with two Lorentzians as one would 
expect [321 06]. The resulting relaxation 
time is shown as mean relaxation time 
(r) = r-r(/3 _1 ) ■ where T is the Gamma 
function, reducing the dependence of r on 
P [13 [29]. 

4. Results 

4-1. Multiple scattering 

This section is organized as follows: After 
looking qualitatively at the data and their 
distortion caused by multiple scattering, 
the two most common ways for the 
quantitative evaluation of QENS data will 
be presented: fits of a sum of Lorentzians 
to the scattering function S^(Q,u) and 
the fit of a stretched exponential to the 
intermediate scattering function I^(Q,t). 

4-1.1. The angular dependence of the scat- 
tered intensity. The algorithm of Paal- 
man and Pings jU] was used to calculate 
the path lengths of the neutrons through 



the samples clS db function of the scatter- 
ing angle. Using the average scattering 
cross section, the probabilities for them 
to be scattered was determined from this 
length [37]. Only those neutrons which 
were scattered exactly once were consid- 
ered to be detected. To facilitate compa- 
rability with the measured data, the re- 
sulting differential scattering cross sections 
were scaled with a factor that was chosen 
such that the calculated integral scattering 
cross section in the covered angular range 
is equal to the measured one for the 0.1 mm 
sample, cf. figured- 

The energy-integrated scattering cross 
section, the diffraction pattern one would 
measure at a standard diffraction instru- 
ment, is shown in figure [3b. Several dif- 
ferences to the calculated one are visi- 
ble: (a) The scattered amount of neutrons 
is much higher than expected from the 
calculation because double scattered neu- 
trons are actually not lost, (b) Increas- 
ing the thickness of the water layer gener- 
ally also increases the scattered intensity. 
This increase is however very small for the 
thick samples - the 3.0 mm sample scat- 
ters hardly more, at small scattering angles 
even less, than the 2.0 mm sample, (c) The 
angular variation of the scattered intensity 
is lower than expected from the calculation. 

The elastically scattered intensity 
determined by integration of an energy 
range of ±60 ueV is shown in figure [3b. 
It can easily be seen that the amount 
of elastically scattered neutrons does not 
scale linearly with the amount of sample 
- the thickest sample even scatters less 
neutrons elastically than the thinnest at 
low angles. 
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Figure 3. From left to right: (a) Calculated angular dependence of the differential 
cross section of H 2 0, caused only by single-scattered neutrons, (b) Measured angular 
dependence of the differential cross section of H 2 0. (c) Measured angular dependence 
of the elastic part of the scattering function of H2O. 
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Figure 4. From left to right: (a) The scattering function S^(Q,uj) as well as the 
imaginary part of the dynamical susceptibility x"{Qj w ) f° r t w0 H2O samples (0.1 mm 
and 3.0 mm sample layer thickness) and values of Q from 0.4 A -1 to 10.0 A -1 in steps 
of 0.2 A -1 . Due to the dynamical range of the spectrometer, the data at low Q do 
not reach high energy transfers and vice versa, (b) The ratio of the normalized double 
differential cross sections of the 3.0 mm measurement and the 0.1mm measurement, 
(c) Shows the same as (b) after conversion from scattering angle 29 to Q. 



4-1.2. Intensity redistribution in (Q,u)- 
space. Comparing the scattering function 
S^(Q,cu) of the thinnest and the thickest 
sample in figure Hk, it becomes clear 
that the elastic signal of the 3.0 mm 
measurement depends much less on Q 
than the one of the 0.1mm measurement, 
in agreement with figure [3b. At higher 
energy transfers, the 3.0 mm sample has 
up to one order of magnitude more signal 
than the 0.1mm measurement. This 
artificial gain in intensity distorts the 
dynamical susceptibility even more than 
the scattering function. 

To look closer at the redistribution of 



the scattered intensity in the quasielastic 
region, the double differential scattering 
cross section of the 3.0 mm measurement 
and the one of the 0.1mm measurement 
were independently normalized to the 
same total scattering cross section in the 
dynamical range. The data of the 3.0 mm 
measurement were then divided by the ones 
of the 0.1mm measurement. The result 
is shown clS db function of energy transfer 
and scattering angle in figure Hb as well 
as a function of energy transfer and Q 
in figure Hb. Values smaller than 10° 
designate where multiple scattering 

caused a loss of intensity, values bigger 
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than 10° those where additional, multiply 
scattered neutrons were detected. 

The pattern shows a slight asymmetry 
with respect to the energy transfer. The 
lower intensity observed on the neutron 
energy loss side was attributed to the 
higher absorption probability for slower 
neutrons. 

At low scattering angles, the signal of 
the 0.1mm measurement decreases more 
rapidly SIS db function of energy transfer 
than the signal of the 3.0 mm measurement 
which is more smeared out, cf. figure UK- 
This leads to the high relative intensity 
at low scattering angles and large energy 
transfers. With increasing scattering angle, 
the shapes of the scattering functions 
become more similar to each other. 

4-1.3. Evaluation of the scattering function 
with two Lorentzians. In the fit with two 
Lorentzians ( [12]) . the properties of the nar- 
row and the broad one, corresponding to 
long-range and localized motions, respec- 
tively, are discussed in the following. 

The arguably most important obser- 
vation in the presented measurements is 
shown in figure Ek and figure Eb: the line 
width of the narrow Lorentzian is largely 
unaffected by multiple scattering. This 
means that a reliable determination of the 
diffusion coefficient is possible even when 
a massive amount of multiple scattering is 
present. This will be detailed in the dis- 
cussion. 

The intensity distribution between 
the two Lorentzians, the EISF, follows 
for the thin samples quite well the 
expected Gaussian-like decrease with Q, 
see figure Eb. However, it is already 
visible in the 0.1mm measurement that 
the EISF does not extrapolate to 1 as Q 



approaches 0. This well-known sign of 
multiple scattering [36J gets much stronger 
with increasing sample layer thickness. 

The width of the broad Lorentzian is 
shown in figure EH. An overall trend to 
larger widths with increasing Q is visible 
on top of which a thicker sample size 
increases the width rather independently 
of Q. As can be seen in figure EK the 
data are at large energy transfers even 
flatter than the fit. This means that 
multiple scattering causes a very broad, 
background-like contribution. 

4-1. 4- Evaluation of the intermediate scat- 
tering function with a stretched exponen- 
tial. The intermediate scattering function 
can be fitted well with a stretched exponen- 
tial as can be seen in figure EK 

The pref actor, displayed in figure Eb, 
is basically the value I^(Q,t = 0), i.e. 
the zeroth Fourier component, which is 
the area under the scattering function, 
jlTmeY S inc(Q^) duJ . The limits of the 
integration are given by the region where 
the scattering function was evaluated - 
which in turn is limited by the dynamical 
range of the experiment. 

The mean relaxation time is influenced 
by multiple scattering, both in its aver- 
age value and its Q dependence. The 
strongest effect can be observed in figure Efc 
at intermediate Q (around 0.7 A -1 ) where 
the apparent speed-up is most pronounced, 
in agreement with Monte Carlo simula- 
tions [TT] . 

The stretching factor is more influ- 
enced than the mean relaxation time |14j . 
Stretching increases (8 decreases) system- 
atically, most at small values of Q, cf. fig- 
ure EH. 
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Figure 5. From left to right: (a) Scattering function of H2O, shown at two Q values 
and two sample layer thicknesses. Also shown is a vanadium measurement as a measure 
of the instrumental resolution and fits of the sum of two Lorentzians to the data (drawn- 
out lines). For the two measurements at low Q, also the two components of the fit 
are shown (dashed lines), (b) The line width of the narrow of the two Lorentzians. 
(c) The elastic incoherent structure factor Aq . (d) The line width of the broad of the 
two Lorentzians. 




Figure 6. From left to right: (a) Intermediate scattering function of H2O, shown at 
two Q values and two sample layer thicknesses. Also shown is the fit of stretched 
exponentials to the data (drawn-out lines), (b) The prefactor a. (c) The mean 
relaxation time multiplied by Q 2 . (d) The stretching exponent /3. 



4-2. Fitting artifacts 

The width of the narrow component of the 
fit to the phospholipid data with a binning 
of SE = lOueV is shown in figure [Q The 
analogous fit to the data with 8E = 2 ueV 
yields the line widths in figure [7_K Also 
shown are in both figures fits of straight 
lines to the medium Q range data points 
and the width of the bins. 

It can be seen that in the case of the 
coarse energy binning shown in figure HJ 
the data points deviate strongly from the 
linear behaviour at low Q and level off to a 
value of about 5E/2. With a finer energy 
binning, figure Efet, they display a much 
more linear behaviour. 

This is shown in detail in figure [7b 



where the fit results for the two energy 
binnings are compared. The effect is 
striking: the levelling off which was 
discussed to be either confinement or a fit 
artifact [TJ is clearly the latter. 

5. Discussion 

5.1. Evaluation of the long-range motion 

Both, geometry and rate of the long- 
range motion are surprisingly little affected 
by multiple scattering. The rate of the 
diffusion is connected to the overall value 
of r narrow /<5 2 or (r) • Q 2 , the geometry to 
the Q-dependence of these quantities. 

The evaluation of the scattering 
function S(Q,u>) with the sum of two 




Figure 7. From left to right: (a) Width of the narrow Gaussian line 
of a three-component fit to the spectra of the fully hydrated phospholipid 
dimyristoylphosphatidylcholine (DMPC). The levelling off of the widths at low 
temperatures and small Q (cf. figure [1]) is practically inexistent. The distance between 
two data points, 5E — 2 ueV, is indicated in the figure. The energy resolution is about 
60 ueV. The solid lines are fits to the filled points in the medium-Q region, (b) Close-up 
view of the low-Q region of the width of the narrow Gaussian line of a three-component 
fit to the spectra of DMPC, measured with different point distances. It can be clearly 
seen that the effect of levelling off at small Q can be nearly completely removed by 
choosing a higher point density. 



Lorentzians where the narrow accounts for 
the long-range motion yields very robust 
parameters for the narrow component. As 
can be seen in figure [5b, only for small Q, 
an increasing scattering power of the sam- 
ple correlates with an increase of the line 
width. 

Of the large deviation of the measured 
line width from the expected values for 
low Q shown in figure [fl up to 125%, 
only a very small part can be attributed 
to multiple scattering: Extrapolating the 
72% and 85% transmission measurements 
to one with 100% transmission, we expect 
a deviation at small Q of about 3%. This 
extrapolation is known to be not suitable 
for the precise evaluation of multiple 
scattering |IOj E] but it shows that the 
effects caused by multiple scattering are 
more than one order of magnitude lower 
than the observed deviations in figure [TJ 



Much in contrast, the fitting artifacts 
distort the observed geometry of the long- 
range motion heavily. It will be shown 
in the following that a systematically 
too large line broadening (or too small 
correlation time) is obtained. 

A Lorentzian line is practically indis- 
tinguishable from a power law for any hw > 
T, namely (T/u) 2 , as shown in figure EJ 
This means that for any energy transfer 
bigger than T, the fit can be shifted to big- 
ger y- values either by increasing the always 
present prefactor a or - completely equiv- 
alent ly - by increasing T. 

Increasing a or increasing T makes a 
difference only at energy transfers smaller 
than T: increasing a increases the y- 
values in this range whereas increasing T 
decreases them. It is therefore this range 
which determines a and T. 

The discrete convolution of the theo- 
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Figure 8. From left to right: (a) Lorentzian curves with different widths plotted 
double-logarithmically. After a plateau, the curves become a power law for hui > L. 
(b) The Riemann sum over a Lorentzian as a function of its width L, shown for different 
values of their centers c. Both, width and center are given in units of the distance 
between two points at which the Lorentzians are evaluated SE. For T < SE/2, the 
Riemann sum deviates significantly from the correct value of 1. The direction of the 
deviation depends on the relative position of the center of the Lorentzian with respect 
to the points. 



retical fit function S^°(Q,uj) with the in- 
strumental resolution R(Q, u) is the origin 
of this problem: whereas R is in the present 
case very well-behaved, S^° becomes so 
sharp peaked as T decreases that the eval- 
uated values at the points near zero en- 
ergy transfer are not representative for the 
whole bin any more. This is because the 
curvature in the center of the Lorentzian is 
— 2/(7iT 3 ) and diverges therefore for r — > 0. 

The consequence is that the integral 
convolution in fl3]) cannot be substituted 
by the discrete convolution fl3j) any more: 
Whereas the integral under the Lorentzian 
is T-independently unity, the Riemann sum 
deviates very strongly when the line width 
is smaller than about half of the distance 
between two evaluated points, T < SE/2. 

It will be shown in the following that 
the systematically too large fit value for 
the line width does not depend on the 
direction of this deviation in the critical 
cases discussed here (T m 1 ueV, SE = 



lOpeV): 

(a) If the sum is calculated too big 
(the center of the Lorentzian is close to 
a measured point), the central points are 
mistakenly calculated too high. Only the 
central point which is evaluated at the 
center of the Lorentzian is in the plateau 
region of the Lorentzian. 

(a.l) The fit can increase T and hereby 
lower the central points by both, the 
natural decrease (figure [8^) and the better 
calculation (figure |8b). Simultaneously, 
the function values in the wings of the 
Lorentzian will increase (figure Ek) which 
will be compensated by a smaller value 
for the prefactor a, making only a small 
counter-contribution to lower the central 
points. 

(a. 2) The fit could alternatively de- 
crease the prefactor to get the central 
points down. This leads to a fit curve that 
is systematically below the data points at 
large energy transfers - it will be shifted 
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upwards by increasing T. 

(b) If the sum is calculated too small 
(the center of the Lorentzian is close to the 
middle in between two measured points), 
the central points are mistakenly calculated 
too low. The points which are evaluated 
close to the center of the Lorentzian are 
already in the power law wing region of the 
Lorentzian. 

(b.l) The fit can increase V and hereby 
increase the central points by both, the 
natural increase (figure [SJi) and the better 
calculation (figure [8b). Simultaneously, 
the function values in the wings of 
the Lorentzian will increase (figure Et) 
which will be compensated by a smaller 
value for the prefactor a which cannot 
compensate the increase caused by the 
better calculation though. 

(b.2) The fit could alternatively in- 
crease the prefactor to get the central 
points up. This leads to a fit curve that 
is systematically above the data points at 
large energy transfers - it would need to be 
shifted downwards by decreasing T. How- 
ever, this would yield even lower central 
points (figure [Hb) and is therefore not an 
option for the fit. 

In conclusion, the fit will always 
yield a line width T that is too broad 
if r > 5E/2. In the example of 
the phospholipid measurements shown in 
figure [1] and figure [3, the deviation from 
a linear behaviour at small Q is decreased 
from 125% (in figure [1]) to about 30% (in 
figure 0) by choosing a smaller energy step 
between measured points. 

This 30% deviation might very well be 
due to the background in the resolution 
measurement, caused by the instrument 
and phonons in the vanadium. Comparing 
the fits using a convolution with the 



instrumental resolution as determined by 
a vanadium measurement at temperatures 
between 4 K and 300 K as well as Gaussians 
with flat background, the fitted values 
differed in this order of magnitude without 
showing clear trends. This will be 
studied in more detail and presented in a 
forthcoming publication. 

5.2. Evaluation of the localized motions 

The evaluation of the localized motions 
in the two-component fit of the scatter- 
ing function is influenced in both, rate 
and geometry. This is because the mul- 
tiply scattered neutrons obtain relatively 
large energy transfers, yielding a nearly 
flat background in the energy direction. 
With an increasing amount of multiple 
scattering, this background increases. The 
broad Lorentzian which accounts for local- 
ized motions therefore broadens and gets 
more intensity. 

The geometry is increasingly hard to 
evaluate from the EISF as it does not 
approach 1 for small momentum transfers 
any more. It seems to be the most 
sensitive indicator for multiple scattering. 
Including an arbitrary scaling factor in the 
definition of the EISF as done previously [1] 
makes an approximate evaluation possible 
at reasonable transmissions [T7] . 

In the following, we will try to explain 
the effects of multiple scattering with 
simple phase space arguments, assuming 
for the calculation of Q that the neutron 
does not change its energy in a scattering 
process and neglecting the geometry of 
the sample [IE]. Knowing that the 
spectrum of multiply scattered neutrons 
is essentially a convolution of the double 
differential cross section with itself for 
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every scattering process jT_2] and that 
two times scattering is the most probable 
multiple scattering event jl9], the Q 
values of the two scattering events are the 
interesting quantities. The width of the 
observed scattering function will then be 
approximatively the sum of the two single 
scattering processes. 

Using Monte Carlo simulations [17], it 
was determined that the probability distri- 
bution of the scattering angles of the sin- 
gle scattering processes is - independently 
of the total scattering angle - nearly trian- 
gular with a maximum probability of 90°, 
corresponding to a Q of already 0.71-Q max 
(elastic), in this case to 1.5 A -1 . A sim- 
ilar assessment has also been used else- 
where [50]. Half of the neutrons will be 
scattered with a larger Q and result in a 
very broad line contributed by the multi- 
ply scattered neutrons. 

The upper limit for the values of Q 
that can possibly participate in a double 
scattering event is given by the wavelength 
of the neutrons: Q™f x = The 
lower limit for Q 12 seen from a given Q tot 
can be deduced from figure where the 
vectorial momentum transfer of the first 
scattering process is denoted with Qx, the 
one of the second with Q 2 , and the total 
(observed) by Q tot — Qi + Q 2 - This 
means for the absolute values: Qi + Q 2 ^ 
Qtot- The distribution of possible Q values 
in the scattering process and therefore 
of relaxation times is therefore largest at 
small Q to t (cf. figure EH). 

The shift with regards to the scat- 
tering angle and the cut-off of Q by the 
choice of the wavelength of the incoming 
neutrons brings about that a change of 
this wavelength influences multiple scatter- 
ing not only for coherently scattering sam- 



ples pUJ EIEJ [TEl EE] but also for incoher- 
ently scattering samples jMllTT], although 
to a smaller extent. 

To understand the resulting line 
width, one needs the conditional probabil- 
ity for the two scattering events given some 
Qtot- We have calculated the mean Q 2 af- 
ter a given Qi as follows (cf. figure Ek): 
The probability for a certain (Qi,Q 2 )- 
combination is directly proportional to the 
number of allowed (Q%, Q 2 )-combinations, 
weighed with the corresponding solid an- 
gles. 

The probability to have an event with 
a certain Q\ )2 is inversely proportional 
to the surface of a sphere with radius 
Qi j2 , 1/(4ttQ1 2 ). Multiplying these two 
probabilities gives the one to reach a 
certain Q tot via a given connection point 
of Qi and Q 2 . Each of the two scattering 
events have to be multiplied with the solid 
angle accessible with a given scattering 
angle, sin(2#i )2 )/A. 

As all connection points which satisfy 
Qi,2 < Qmax can be used equivalently, this 
probability has to be multiplied with the 
circumference / of the circle around Qtot, 

i = 7^-Mi-(Qtot-Qi-Q2)- 

Qtot *> 

(Qtot - Qi + Q 2 ) • 
(Qtot + Qi - Q 2 ) ■ 
(Qtot + Qi + Q 2 )] 1/2 } , (14) 

where 5i denotes the real part of the 
argument. 

The result, symmetrical in Qi and Q 2 , 
is displayed in figure [9b. It can be seen that 
even if one scattering process has already 
Qi = Qtot, the other scattering process still 
has a nonzero Q 2 . Two limiting behaviours 
can be observed: On the one hand, for 
Qi and Q tot > 1.0 A -1 , the mean Q 2 
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Figure 9. From left to right: (a) Two scattering processes, denoted with subscript 1 
and 2, with their respective scattering vectors Qi and Q 2 combine to a total scattering 
process with Q to t • The black dash-dotted lines denote the external borders in between 
which a connection point between the two scattering events must be so that the given 
Qtot can be reached with two quasi-elastic scattering processes. Red dashed line: 
section of the sphere on which Q^ may end and the paper plane. Note that only points 
on this sphere which are within the black dash-dotted line can contribute to the shown 
scattering process. Blue dotted line: projection of the circle (which is perpendicular 
to the paper plane) on which the connection point between the two Q" can revolve, 
(b) Mean Q 2 for a given combination of Qi and Qtot- 



depends basically only on |Q tot — Qi\ and 
is never less than about 1.0 A -1 . On the 
other hand, for Qi and Qtot ^$ 1.0 A -1 , we 
observe Q 2 ~ max (Qi, Qtot)- Scattering 
events where the second scattering process 
has a very small Q compared to Q\ 
and Q t ot, i.e. broadening the signal only 
slightly, are very rare. This causes the 
distinctness between the signal and the 
much broader contribution of multiply 
scattered neutrons. 

5.3. Phospholipid motions 

It is clear that the long-range motion ex- 
tracted from a measurement with an obser- 
vation time of 60 ps is not the macroscopic 
long-range motion. This manifests for ex- 
ample in a higher velocity [T] and lower 
influence of additives [51 J than observed 
macroscopically. Although the molecule 



travelled only a distance of about 1 A [TJ 
during this short time, the interactions 
with the neighbours, the cage, are already 
established. The observed flow-like mo- 
tion is not the initial phase of free flight 
before the molecule experiences the pres- 
ence of the surrounding neighbours. This 
is evidenced by both, the clear existence 
of localized (caged) motions and the small 
flow velocity which is nearly three orders of 
magnitude below the thermal expectation 
value pp. 

As shown in the present contribution, 
the neutron scattering data show no 
significant evidence for a restriction of 
the flow-like long-range motion to a cage 
within the observation time - above and 
below the main phase transition. 

We explain this with "floppy modes": 
the molecules are pushed around by density 
fluctuations resembling longitudinal acous- 
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tic phonons in crystals. Below the phase 
transition, a restoring force will eventu- 
ally push them back to the original place. 
Due to the "Soppiness" of the material, the 
molecules travel more than the 1 A dur- 
ing which the experiment follows them be- 
fore the restoring force gets big enough to 
restrict their motion. Therefore, this ef- 
fect does not leave a trace in the data. 
Above the phase transition, the molecules 
are transported by many short subsequent 
uncorrelated flow events. 

These low-frequency "floppy modes" 
are known from disordered colloidal crys- 
tals |52j, colloidal glass formers, and su- 
percooled liquids and signal the onset of 
macroscopic elasticity [53J. Phospholipid 
membranes also exhibit elasticity above 
and below the main phase transition [31] 
and the observed flow patterns [2] resemble 
very much dynamical heterogeneities [STj, 
a hallmark of glassy dynamics, which were 
again suggested to be linked with the 
"floppy modes" [55]. 

The reason that the dynamics of phos- 
pholipids resembles so much dense, super- 
cooled systems is probably the very high 
density in the head group area of the phos- 
pholipids [5T]- This similarity has proven 
to be extremely useful already in the past 
when the free volume theory developed for 
glasses [52] was adapted to membranes [57J, 
135] . Transposing the newer concepts de- 
veloped in glass physics, dynamical het- 
erogeneities and "floppy modes", to phos- 
pholipid membranes seems not only con- 
sequent but can also explain the observed 
features effortless. 



6. Summary and conclusions 

The molecular motion of the phospholipid 
DMPC was studied below and above the 
main phase transition on a time scale 
of about 60 ps, taking the influence of 
multiple scattering and fit artifacts into 
particular consideration. 

Multiple scattering was shown to be 
only a minor nuisance for the determina- 
tion of the long-range motion of an inco- 
herent scatterer. Exemplary, H2O samples 
with transmissions between 85% and as low 
as 2% were studied. The multiply scattered 
neutrons are lost in the elastic region of the 
signal and create a background which is a 
very broad line in energy space. In a mul- 
ticomponent analysis, they affect therefore 
only the evaluation of localized motions. 
This, together with the excellent agreement 
of the results with Monte Carlo simulations 
shows that one could more often dare to use 
strongly scattering samples for quasielastic 
measurements. 

The fit of a theoretical scattering 
function to measured data after numerical 
convolution with a measured resolution 
function can give erroneous results. A 
delicate combination of correlated fitting 
parameters and a discretization problem 
during the convolution with the measured 
instrumental resolution causes an error 
in an unexpected place: the line width. 
Analyzing the data with different bin 
widths should give indications if this 
problem is present. The lower limit for 
the bin width is often given by the online 
histogramming of the neutrons during the 
experiment - a complete event recording of 
the detected neutrons, in principle possible, 
would remove this limitation. 

Flow-like motions of the phospholipid 



Slow short-time phospholipid motions, multiple scattering, and fitting artifacts 17 



molecules bear a striking resemblance to 
dynamical heterogeneities above the main 
phase transition and to "floppy modes" 
below. The confinement of the molecules 
below the main phase transition is not 
visible in the long-range component of the 
fit due to the short observation time of the 
experiment. It seems that the description 
of motions in the phospholipid membrane 
could much profit from transferring the up- 
to-date concepts used in glass physics. 
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